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Abstract 



The search for D°D° mixing carries a large discovery potential for new physics since the 
D°D° mixing rate is expected to be very small in the Standard Model. The past decade has seen 
significant experimental progress in sensitivity, from 20% down to 0.37%. This paper discusses 
the techniques, current experimental status, and future prospects for the mixing search. Some 
new ideas, applicable to future mixing searches, are introduced. The conclusion is that while it 
is possible that the mixing sensitivity may decrease to 10~ 5 around the year 2000, reaching the 
' 10~ 6 level will certainly be quite difficult. 
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1 Introduction 



Following the discovery of the D° meson at SPEAR in 1976, experimenters began to search for D°D° 
mixing, using a variety of techniques. The past decade has seen significant experimental progress in sensitivity 
' (from 20% to 0.37% Q to p3[), as can be seen in Figure 1. Much of the enthusiasm for searching for D°D° 

mixing stems from the belief that the search carries a large discovery potential for New Physics, since the 
mixing rate R m i x in g = B(D° — » D° — > /)/ B(D° — > /) is expected to be very small in the Standard Model. 
One can characterize D°D° mixing in terms of two dimensionless variables: x = Sm/j^ and y — 7-/7+, 
where the quantities j± and 8m are defined by j± = (71 ± 72)/2 and dm = mi — m\ with m^, 7, [i = 1,2) 
being the masses and decay rates of the two CP (even and odd) eigenstates. Assuming a small mixing, 
namely, 5m, -C 7+ or x,y <C 1, we have R m ixing = {x 2 + 2/ 2 )/2. Mixing can be caused either by x 7^ 
(meaning that mixing is genuinely caused by the D° — D° transition) or by y 7^ (meaning mixing is caused by 
the fact that the fast decaying component quickly disapears, leaving the slow decaying component which is a 
mixture of D° and D°). Theoretical calculations of D°D° mixing in the Standard Model are plagued by large 
uncertainties. While short distance effects from box diagrams are known |l4| to give a negligible contribution 
(~ 10 -10 ), the long distance effects from second-order weak interactions with mesonic intermediate states 
may give a much larger contribution. Estimates of R m i x ing from long distance effects range from 10 -7 to 
10~ 3 fig] . However, it has recently been argued by Georgi and others that the long distance contributions are 
smaller than previously estimated, implying that cancellations occur between contributions from different 
classes of intermediate mesonic states jL(|, and the prevailing conclusion within the Standard Model seems 
to be that R m i X ing < 10~ 7 jl7|. A measurement of such a small mixing rate is not possible with present 
experimental sensitivity. However, the observation of a larger value for R m i x ing caused by x 7^ would imply 
the existence of new physics beyond the Standard Model [l8| . Examples includes flavor-changing neutral 
currents mediated by the exchange of a non-standard Higgs scalar with a mass of a few TeV/c 2 , which could 
lead to Rmixing as large as 0.5%. 
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Recently, CLEO has observed a signal for D° -> K+ir- @, and found R = B{D° -> K+tt')/ B(D° -> 
K-TT+) ~ 0.8% @. Normally, D a decays by Cabibbo favored decay D° -> #"71-+ and L>° -> i^+Tr^. A 
signal for D° — > X + 7r~ could indicate mixing of D° — > D°. But it could also indicate a different decay 
channel, namely, Doubly Cabibbo Suppressed Decay(DCSD) D° — > if + 7r~, which is suppressed with respect 
to the Cabibbo favored decay by a factor of tan^Oc ~ 0.3% where 0c is the Cabibbo angle. Unfortunately, 
without a precision vertex detector, CLEO is unable to distinguish a potential mixing signal from DCSD. As 



Purohit pointed out |21 : "The CLEO II result makes the entire subject of D°D° mixing very interesting. 
It really calls for a fixed-target experiment to use its decay time resolution to decide whether the signal is 
due to DCSD or mixing" . If the number of reconstructed charm decays can reach 10 s around the year 2000, 
that would allow one to reach a new threshold of sensitivity to D°D° mixing, and perhaps actually observe 
it. This is why charm mixing has been singled out for its own working group at this workshop. 

This paper is organized as follows: Section 2 discusses the techniques which can be used to search 
for mixing, including two new ideas. One of them is to use D° — ► K + K~ , 7r + 7r~ etc. to study mixing (see 
2.1.1), and the other is to use the difference in the resonant substruture in D° — > K + ir~ir°,D — > K + tt~ ir + ir~ 
etc. to distingusih mixing and DCSD (see 2.1.2). In each case, the relevant phenomenology will be briefly 
presented. Section 3 discusses the present status and future prospects of searching for mixing at different 
experiments. In section 4, a comparison of the future prospects of the different experiments with different 
techniques will be given. A brief summary is in Section 5. 



2 The Techniques 



The techniques which can be used to search for mixing can be roughly divided into two classes: hadronic 
and semi-leptonic. Each method has advantages and limitations, which are described below. 



2.1 Hadronic method 

The hadronic method is to search for the D° decays D° — > K + ir~ (X). These decays can occur either 
through D°D° mixing followed by Cabibbo favored decay D° -> D° -> K+tt-(X), or through DCSD 
D° — > K + ir~ (X). This means that the major complication for this method is the need to distinguish 
between DCSD and mixing ]22]| . The hadronic method can therefore be classified according to how DCSD 
and mixing are distinguished. In principle, there are at least three different ways to distinguish between 
DCSD and mixing candidates experimentally: (A) use the difference in the decay time-dependence; (B) use 
the possible difference in the resonant substructure in D° — ► K + ir~ it® ,K + ir~ n + ir~ : etc. modes; (C) use the 
quantum statistics of the production and the decay processes. 

Method (A) requires that the D° be highly boosted and so that the decay time information can be 
measured. Method (B) requires knowledge of the resonant substructure of the DCSD decays, which is 
unfortunately something about which we have no idea at this time. Finally, method (C) requires that one 
use e + e~ annihilation in the charm threshold region. In the following, we will discuss these three methods 
in some detail. 



2.1.1 Method A -use the difference in the time-dependence of the decay 

This method [^3) is to measure the decay time of the D° — > K + ir~ decay. Here the D° tagging is 
usually done by using the decay chain D* + — > D°tt+ followed by D° — > K + ir~ . The 7r+ from D* + has a soft 
momentum spectrum and is refered to as the slow pion. The idea is to search for the wrong sign D* + decays, 
where the slow pion has the same charge as the kaon arising from the D° decay. This technique utilizes the 
following facts: (1) DCSD and mixing have different decay time-dependence, which will be described below. 

(2) The charge of the slow pion is correlated with the charm quantum number of the D° meson and thus 
can be used to tag whether a D° or D° meson was produced in the decay D* + — > D°ir+ or D*~ — * D°tt~ . 

(3) The small Q value of the D* + decay results in a very good mass resolution in the mass difference 
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AM = M(D* + ) - M(D°) - M(ir+) and allows a D* + signal to obtained with very low background. (4) 
The right sign signal D* + — > -D°7r+ followed by D° — > K~n + can be used to provide a model- independent 
normalization for the mixing measurement. 

A pure D° state generated at t — decays to the K + tt~ state either by D°D° mixing or by DCSD, 
and the two amplitudes may interfere. The amplitude for a D° decays to K + tt~ relative to the amplitude 
for a D° decays to K~tt + is given by 



A = ^R m bdng/2 t + v/Rdcsd e*+ (1) 

where <p is an unknown phase, t is measured in units of average D° lifetime. Here Rdcsd = \p\ 2 where p is 
defined as p = Amp(D° — > K + Tr~)/Amp(D° — > i"T + 7r~), denoting the relative strength of DCSD. We have 
also assumed a small mixing; namely, 5m, 7- <C 7+ or x, y -C 1, and CP conservation. 

The first term, which is proportional to t, is due to mixing and the second term is due to DCSD. It 
is this unique attribute of the decay time-dependence of mixing which can be used to distinguish between 
DCSD and mixing. Now we have: 



l(D° K+ir-)(t) OC (RdCSD + v^RmixingRDCSD tC0S(j) + 1/2 Rmbdng^K* (2) 

Define a = R m ixing/RDCSD , which describes the strength of mixing relative to DCSD. Equation (2) can then 
be rewritten as: 

K7) -► K+n-)(t) oc R DC sd(1 + v 7 ^ tcoscj) + l/2at 2 )e~* (3) 

From this equation, one may read off the following properties: (1) The mixing term peaks at t = 2. 
(2) The interference term peaks at t = 1. (3) A small mixing signature can be greatly enhanced by DCSD 
through interference (with coscp ^ 0) at lower decay times. The ratio between the interference term and 
the mixing term, denoted £(t), is given by £(t) = y/8/a coscp/t oc y/l/a. So when a — > 0, £ — > oo. 
For example, with coscj)— I, at t = 1 for a = 10%, 1%, 0.1%, 0.01%, 0.001% (corresponding to R m ixin g = 
10- 3 ,10~ 4 ,10- 5 ,10- 6 ,10- 7 ) one has = 9,28,90,280,900 respectively. (4) Only for t > ^/ff^\cos(j)\ 
does the interference term become smaller than the mixing term. (5) I(to) = happens and only happens 
when coscj) = —1, and only at location to = \j2ja. For a = 10%, 1%, 0.1%, one has to = 4.5,14.1,44.7. 
(6) One can obtain a very pure DCSD sample by cutting at low decay time: t < ( ~ 0.1. At such low t, 
the mixing term drops out and leaving only the interference term. Let's define the purity of DCSD to be 
Vdcsd = 1 - Jo V^a coscPite-*)/ J Q C e'K For C ~ 0.1 one get Vdcsd = V^acoscj) f/(2 + C)- Let's take 
C = 0.2, for a = 10%, 1%, we can get Vdcsd = 96%, 99% pure DCSD respectively. (7) The cut t < ( cuts 
off only 1 — (1 + C) e ^ ~ C 2 /2 = 2% from the whole interference term. 

While Property (1) tells us that the mixing term does live at longer decay time, Property (3) tells us 
clearly that we should not ignore the interference term. In fact, that's the last thing one wants to ignore! 
(unless we know for sure cos cf> = 0). It is the commonly believed "annoying background", namely DCSD, 
that could greatly enhance the chance of seeing a very small mixing signal. In semi-leptonic method, one 
does not have this advantage. For a very small mixing rate, almost all the mixing signature could show up 
in the interference term, not in the mixing term, as long as cos cp ^ 0. Property (2) tells us at which location 
one expect to find the richest signature of a potential small mixing, which is where the interference term 
peaks: t ~ 1 (why should one keep worrying about long lived DCSD tails? let's hope for cos0 ^ first.) 
Property (5) shows that destructive interference is not necessarily a bad thing. In fact, it could provide extra 
information. For example, if cos cp = — 1, then one should find I(to) = at to = \/2/a, see Figure 5. For the 
general case, interference will lead to very characteristic time distribution, as can be clearly seen in Figure 
6. Properties (6) and (7) show that we can study DCSD well without being confused by the possible mixing 
component. This will become important when we discuss method B. 

Therefore the signature of mixing is a deviation from a perfect exponential time distribution with the 
slope of 7 + f]. Our ability to observe this signature depends on the number of D° — > K + n~ events we will 

J One can use D° — > K~n + to study the acceptance function versus decay time. 
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have. Right now this is limited by the rather poor statistics. Figures 3 and 4 show each term with a = 10% 
and cos<p = ±1 (with Rdcsd = 1)- 

It is interesting to point out here that there is also a possibility, previously unrecognized, of using the 
Singly Cabibbo Suppressed Decays (SCSD), such as D° — > K + K~, tt + it~ to study mixing. This is because 
(assuming CP conservation) those decays occur only through the CP even eigenstate, which means the decay 
time distribution is a perfect exponential with the slope of 71. Therefore, one can use those modes to measure 
71. The mixing signature is not a deviation from a perfect exponential (again assuming CP conservation), 
but rather a deviation of the slope from (71 + 72V2. Since 7+ = (71+72)72 can be measured by using the 
D° —y K~n + decay time distribution, one can then derive y — 7-/7+ = (71 — 72)/ (71 +72)- Observation 
of a non-zero y would demonstrate mixing caused by the decay rate difference (Rmixing = [x 2 + y 2 )/2). 
It is worth pointing out that in this case other CP even (odd) final states such as D° — > K + K~ K + K~ 
(K + K~ir°) can be also used to measure 71(72)- In addition, there is no need to tag the D° nor know the 
primary vertex location, since we only need to determine the slope. I should point out also that this method 
is only sensitive to mixing caused by the decay rate difference between the two eigen states, not to mixing 
caused by the mass difference x = Sm/^ + (5m — m 2 — mi). The sensitivity of this method is discussed in 
section 4.1. 



2.1.2 Method B - use difference in resonance substructure 

The idea of this new method is to use the wrong sign decay D* + — > D w+ followed by D° — » K + ir~ir°, 
K + 7r~7r + 7r~, etc., and use the possible differences of the resonant substructure between mixing and DCSD 
to study mixing. There are good reasons to believe that the resonant substructure of DCSD decay is different 
from that of mixing (Cabibbo favored decay, CFD). We can use the D° — > K + ir~ir° decay as an example. 
For CFD and DCSD, the true yield density n(p) at a point p in the Dalitz plot can be written as: 

n[p) cx |/x e^Asb + f 2 e^BW p+ (p) + / 3 e**BW K .- (p) + h e»*BW R . ( P )\ 2 (4) 

where /, are the relative amplitudes for each component and <f>i are the interference phases between each 
submode. is the S-wave three-body decay amplitude, which is flat across the Dalitz plot. The various 
BW(p) terms are Breit-Wigner amplitudes for the K*ir and Kp sub-reactions. Note that in general 0, 

f DCSD /f CFD f DCSD /f CFD ( . j} (5) 

This means that the resonant substructure (the true yield density n(j>)) for DCSD is different from that of 
mixing |. As both DCSD and mixing contribute to the wrong sign decay, the yield density for the wrong 
sign events n w (p) will have a complicated form, due to the fact that for each submode DCSD and mixing 
may interfere with each other. Just like in method A, for very small mixing, the interference term between 
DCSD and mixing could be the most important one. 

In principle, one can use the difference between the resonant substructure for DCSD and mixing events to 
distinguish mixing from DCSD. For instance, combined with method A, one can perform a multi-dimensional 
fit to the data by using the information on AM, M(D°), proper decay time t and the yield density on Dalitz 



It has been pointed out 24, [34| that it is unlikely that just one universal suppression factor will affect 
the individual DCSD. For example, SU(3) breaking can introduce a significant enhancement for D° DCSD 
decays D° — > K + ir~ , K* + ir~ , while SU(6) breaking can introduce a sizeable suppression relative to the naive 
expectation for L>° DCSD decay D° -> K + p~. 

3 The sign of the interference between each submode changes whenever cosOr (Or is the helicity angle of 
the resonance) changes sign. This is the same for both the Cabibbo favored decay and the DCSD. This can 
be easily seen from the Breit-Wigner amplitude which describes the strong resonances and decay angular 
momentum conservation: BWn oc -77 — 1°/ eR ^ — 7^ where Mr and Tr are the mass and width of the Mn 
resonance (the K* or p). The difference between Cabibbo favored decay and DCSD comes from the relative 
amplitude fi for each submode and the interference phase term e l ^ i . 
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plot n w (p, t). The extra information on the resonant substructure will, in principle, put a much better 
constraint on the amount of mixing. Of course, precise knowledge of the resonant substructure for DCSD 
is needed here and so far we do not know anything about it. Because of this, for current experiments this 
method is more likely to be a complication rather than a better method when one tries to apply method 
A to D° — ► K + ir~Tr° (see section 3.2 and pj|) or D° — ► K + ir~ir + Tr~ . In principle, however, one can use 
wrong sign samples at low decay time (which is almost pure DCSD, see section 2.1.1.) to study the resonant 
substructure of the DCSD decays. In the near future, we should have a good understanding of DCSD decays 
and this method could become a feasible way to search for mixing. 



2.1.3 Method C — use quantum statistics of the production and decay processes 

This method is to search for dual identical two-body hadronic decays in e + e~ — > — > D°D°, such as 
(K ~tt + )(K ~7r + ), as was first suggested by Yamamoto in his Ph.D thesis (26). The idea is that when D°D° 
pairs are generated in a state of odd orbital angular momentum (such as ^S"), the DCSD contribution to 
identical two-body pseudo-scalar-vector (D — > PV) and pseudo-scalar-pseudo-scalar (D — > PP) hadronic 
decays (such as (K ~tt + )(K~tt + )) cancels out, leaving only the contribution of mixing [^6[ [24|, |27j. As many 
people have asked about this, I would like to show here the essence of Yamamoto's original calculation for the 
{K-tt+){K-it + )_ case. Let's define e t {t) = e -i"n*-7it/a (j, = 1,2) ande±(f) = {e 1 (t)±e 2 {t))/2. A state that is 
purelyjD ) or \D°) at time t = will evolve to \D(t)) or \D(t)) at time t, with \D(t)) = e + {t)\D°)+e-(t)\D°) 
and_|.D(i)) = e_(t)|D°) + e+(i)|5°). In e+e" -> -> D°Z) , the £>°l) pair is generated in the state 
D Q D° — D°D° as the relative orbital angular momentum of the pair £ = 1. Therefore, the time evolution 
of this state is given by |D(t)5(f)) - |L>(i)L>(^)), where t (f) is the time of decay of the D (5). Now the 
double-time amplitude A w (t,t') that the left side decays to K~ir + at i and the right side decays to K~ir + 
at t' , giving a wrong sign event (K ~n + )(K~Tr + ) 1 is given by: 

A w {t,t') = (e+(t)e_(f ) - e_(i) e+ (i'))(a 2 - 6 2 ) (7) 

where a = (K~tt + \D°) is the amplitude of the Cabibbo favored decay D° — > K~ir + , while 6 = (if _ 7r + \D°) 
is the amplitude of DCSD D° — * K~tt + . Similarly, the double-time amplitude „4 r (i, i') for the right sign 
event (K~ir + )(K + Tr~) is given by: 

A(M') - (e+(t)e+(0 - e_(t) e _(t'))(a 2 - 6 2 ) (8) 

One measures the wrong sign versus right sign ratio R, which is: 

N(K-tt+,K-tt+) +N(K + n~,K+7r-) _ JJ A w (t,t') dt dt' 
~ N{K-ir+,K+ir-) +N(K+Tr-,K-n+) ~ JJ A r (t,t') dtdt 1 ^' 

Note in taking the ratio, the amplitude term (a 2 — b 2 ) in Equations (7) and (8) drop out. Thus, clearly R 
does not depend on whether b is zero (no DCSD) or finite (with DCSD). Integrating over all times, one then 
obtains R = {x 2 + y 2 )/2 = R m i X mg, where x and y are defined as before. 

This is probably the best way to separate DCSD and mixing. The exclusive nature of the produc- 
tion guarantees both low combinatoric backgrounds and production kinematics essential for background 
rejection. This method requires one use e + e~ annihilation in the charm threshold region. Here the 
best final state is (K~tt + )(K~tt + ). In principle, one can also use final states like (K ~ p + )(K~ p + ) or 
(K*~ir + )(K*~Tr + ), etc., although again there are complications. For example, it is hard to differentiate 
experimentally (K~ p + )(K~ p + ) from (K~ p + )(K~ir + Tr°), where DCSD can contribute. With high statistics, 
in principle, this method could be combined with method B. 

It has been pointed out that quantum statistics yield different correlations for the D°D° decays from 
e+e- -> D a D a ,D a D a j,D°D°ir a fjf|. The well-defined coherent quantum states of the D°D° can be, in 
principle, used to provide valuable cross checks on systematic uncertainties and, more importantly, to extract 
x = 8m 1 7 + and y — 7-/7+ (which requires running at different energies) if mixing is observed []33"[ . 
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2.2 Semi-leptonic method 

The semi-leptonic method is to search for D° — ► D° — > Xl~v decays, where there is no DCSD involved. 
However, it usually (not always!) suffers from a large background due to the missing neutrino, in addition, 
the need to understand the large background often introduces model dependence. In the early days, the small 
size of fully reconstructed samples of exclusive D° hadronic decays and the lack of the decay time information 
made it difficult to constrain the D°D° mixing rate using the hadronic method, many experiments used semi- 
leptonic decays. The techniques that were used were similar — searching for like-sign or fj,~fi~ pairs 
in fi+N -> fj+{fi+fj+)X @, § and tT Fe -> fi+fi+ §, ir~W -> p+ fj+ @. These techniques rely on the 
assumptions on production mechanisms, and the accuracy of Monte Carlo simulations to determine the large 
conventional sources of background. 

There are other ways of using the semi-leptonic method. The best place to use the semi-leptonic method 
is probably in e + e~ annihilation near the charm threshold region. The idea is to search for e + e~ — » — » 
D°D° -> {K-l+v){K-l+v) or e+e" -> D~ D*+ -> {K+tt~tt~)(K+1~v)tt+ g The latter is probably 
the only place where the semi-leptonic method does not suffer from a large background. It should have a 
low background, as there is only one neutrino missing in the entire event, threshold kinematics constraints 
should provide clean signal. 

It has been pointed out that one can not claim a Z?°Z?° mixing signal based on the semi-leptonic method 
alone Q (unless with the information on decay time of D°). Nevertheless, one can always use this method to 
set upper limit for mixing. 



3 Mixing Searches at Different Experiments 

3.1 e+e running on *"(3770) MARK III, BES, Tau-charm factory 

The MARK III collaboration was the first (though hopefully not the last) to use the e + e~ — > VP" — * 
D°D° technique. They reported preliminary results for two "wrong-sign" D° decay events (unpublished) [ pT| , 
one was consistent with K~ p + vs. K~ p + , while the other one was consistent with K*°ir° vs. K*°ir°. This 
was a very interesting result at that time, and had a strong influence on the subject. However, one cannot 
draw a firm conclusion about the existence of D°D° mixing based on these events. There are at least two 
reasons: (I) The background study has to rely on Monte Carlo simulation of the PID (particle identification 
- Time-of-Flight) . As Gladding has pointed out: "These results must be considered preliminary because the 
calculation of the confidence level is sensitive to the tails of PID distribution for the background" 63]; (2) 
Assuming that the Monte Carlo background study is correct, and that the events are real, one still cannot 
claim the two events are due to mixing, for example, the non-resonant decays D° — ► Ktttt may contribute 
to one side of the pair in each of the events, in which DCSD can contribute. 

The MARK III puzzle can be completely solved at a r-charm factory, which is a high luminosity 
(10 cm s ) e + e _ storage ring operating at center-of-mass energies in the range 3-5 GeV. The perspectives 
for a D°D° mixing search at a r-charm factory have been studied in some detail |^8|, ^9). I will outline here 
the most important parts. The best way to search for mixing is probably to use e + e~ — * ty" — ► D°D° — ► 
(K~ir + )(K~Tr + ). The sensitivity is not hard to estimate. Assuming a one year run with a luminosity of 
10 33 cm -2 s -1 , 5.8 x I0 7 D°s would be produced from Therefore about 9 x I0 4 (K~tt + )(K + %~) events 
would be produced. About 40% of them (3.6 x 10 4 ) could be fully reconstructed. A detailed study has 
shown that the potential dominant background comes from doubly misidentified (K~ir + )(K + Tr~), and if 
TOF resolution is 120 ps, this background can be kept to the level of one event or less. This means one can 
set an upper limit at the 10 -4 level. 

4 Bigi |}3| has pointed out that an observation of a signal on D° — > l~X establishes only that a cer- 
tain selection rule is violated in processes where the charm quantum number is changed, namely the rule 
ACharm = — AQi where Qi denotes leptonic charge. This violation can occur either through D°D° mixing 
(with the unique attribute of the decay time-dependence of mixing) , or through new physics beyond the 
Standard Model (which could be independent of time). 
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As mentioned section 2.2, the best place to use the semi-leptonic method is probably at a r-charm 
factory. One good example is to search for e + e~ — > D~D* + — > (K + it~it + )(K + 1~i')it+ . It is expected that 
this method can also have a sensitivity at the 10~ 4 level. There are many other independent techniques 
that one can use for a mixing search at a r-charm factory. By combining several independent techniques 
(which require running at different energies), it was claimed that D°D° mixing at the 1CP 5 level could be 
observable p9| . 

There have been several schemes around the world for building a r-charm factory. If such a machine 
is built, it could be a good place to study mixing. At the workshop, Walter Toki told us the history of the 
r-charm factory: one was proposed at SLAC in 1989 and at Spain in 1993. It was discussed at Dubna in 
1991, at IHEP (China), and at Argonne (this workshop). It will be discussed again at IHEP (China) soon. 
Let us hope that we will have one in the not too distant future. 



3.2 e+e running near T(45) ARGUS, CLEOII, CLEO III, B factory 

Using the CLEO II data sample, with an integrated luminosity of 1.8 fb" 1 at and near the T(4S) 
resonance, CLEO has observed a signal for D° — > K + ir~ p(| from the decay chain D* + — > -D°7r+ — * 
(JC + 7r")7r+, as can be seen in Figure 2. 

Without a precision vertex detector, CLEO II can only in effect measure the rate B(D° — > Ktt) 
integrated over all times of a pure D° decaying to a final state Ktt. The ratio K=B(D° — > K + n~) 
/B(D° — > K~ir + ) is given by integrating equation (2) over all times: 



R = R-mixing + R-DCSD + V^R•mixingRDCSDCOS < ?! , (10) 

CLEO II finds R = (0.77 ± 0.25 (stat.) ± 0.25 (sys.))%. This signal could mean one of two things: (1) 
mixing could be quite large, which would imply that mixing can be observed in the near future; (2) the 
signal is dominated by DCSD. The theoretical prediction for Rdcsd is about 2 tan 4 6c ~ 0.6% [^4[ |35| , 
which is quite consistent with the measured value. It is, therefore, believed by many that the signal is due 
to DCSD, although it remains consistent with the current best experimental upper limits on mixing, which 
are (0.37 - 0.7)% g§ and 0.56% Q. 

CLEO has also tried to use hadronic method B, by searching for D° — > K + ir~Tr°. The excellent photon 
detection at CLEO II allows one to study this mode with a sensitivity close to D° — > K + ir~ mode. The main 
complication faced here is that (as discussed in section 2.1.2) the resonant substructure is not necessarily the 
same for wrong sign and right sign decays. Because of this, the interpretation of R as Rmixing or Rdcsd will 
be complicated by the lack of knowledge of the details of the interference between submodes (and also the 
decay time information). Moreover, one has to worry about the detection efficiency across the Dalitz plot. 
Setting an upper limit for each submode is clearly very difficult. CLEO has thus set an upper limit on the 
inclusive rate for D° -> K+tt-tt as R = B(D° -> K+tt-tt ) /B(D° -> R-tt+tt ) < 0.68% @. Note this 
upper limit includes the possible effects of the interference between the DCSD and mixing for each submode 
as well as the interference between submodes. 

This summer, CLEO will install a silicon vertex detector (SVX) with a longitudinal resolution on vertex 
separation around 75 /im. This will enable CLEO to measure the decay time of the D°, and reduce the 
random slow pion background (as the resolution of the D* + - D° mass difference is dominated by the angular 
resolution on the slow pion, this should be greatly improved by the use of the SVX). By the year 2000, with 
CLEO III (a symmetric B factory) and asymmetric B factories at SLAC and KEK, each should have about 
thousands D° -> K + K-(X),n + ir-(X) and a few hundred D° -> K + tt~ (and perhaps D° -> K+tt-tt , 
K + ir~n + ir~ too) signal events with decay time information for one year of running. The typical decay 
length of D° (£) is about a few hundred ^m, and the resolution of the decay length (ac) is about 80 fim 
(£/<7£ ~ 3). The sensitivity to mixing at CLEO III and asymmetric B factory has not been carefully studied 
yet. A reasonable guess is that it could be as low as 10~ 4 . If mixing is indeed as large as DCSD, it should 
be observed by then. 



7 



3.3 Fixed target experiments E615,E691, E791, E687,E831 



A significant amount of our knowledge has been gained from Fermilab fixed target experiments, and in 
fact the current best upper limits on mixing have emerged from these experiments (E615, E691), and will 
come from their successors E687, E791 and E831 soon. 

The best upper limit using the semi-leptonic method comes from the Fermilab experiment E615, which 
used a 255 GeV pion beam on a tungsten target. The technique is to search for the reaction irN — > D°D° — > 
(K~ n + is)D Q — > (K~ fi + v)(K~ fi + u), where only the final state muons are detected (like sign fi + fi + or \i~ yT 
pairs). Assuming a{cc) ~ A 1 nuclear dependence, they obtained R m i x i ng < 0.56% |lC(] . 

The best upper limit using the hadronic method by measuring the decay time information comes from 
E691, which is the first high statistics fixed target (photoproduction) experiment. In fact, E691 was the first 
experiment which used the decay time information (obtained from the excellent decay time resolution of their 
silicon detectors) to distinguish DCSD and mixing. The decay chains D* + — > D°ir^ followed by D° — > K + n~ , 
K + Tr~7r + ir~ were used. The upper limits from the D° — > K + tt~ mode are R m i x ing < (0.5 — 0.9)% and 
Rdcsd < (1.5 - 4.9)% , while the upper limits from D° — > K+tt^tt+tt^ are Rmixing < (0.4 - 0.7)% and 
Rdcsd < (1.8 — 3.3)% . The combined result gives R m ixin g < (0.37 — 0.7)%. The ranges above reflect 
the possible effects of interference between DCSD and mixing with an unknown phase(0). Note that for 
D° — ► K + ir~ir + ir~ , the resonant substructure in the Cabibbo favored and DCSD decays has been ignored. 

At this workshop, both E687 and E791 have reported their preliminary result from part of their data. 
One can find the details in Jim Wiss's talk. The best upper limits on mixing should come from these 
two experiments soon. It is worth pointing out here that both the E687 and E791 results reported in the 
workshop are based on the assumption that there is no interference between DCSD and mixing. Future 
analysis should include the interference term for the reasons discussed in section 2.1.1. 



4 Comparison of Different Experiments 

4.1 Hadronic method A 

This measurement requires: (1) excellent vertexing capabilities, at least good enough to see the inter- 
ference structure; (2) low background around the primary vertex. The background level around the primary 
vertex could be an important issue as the interference term in equation (2) does peak at t = 1. In addi- 
tion, low background around primary vertex means that one does not suffer much from random slow pion 
background and also one can measure the DCSD component at low decay time well. This is important for 
understanding DCSD at large decay times. The vertexing capabilities at e + e~ experiments (C/a ~ 3) for 
CLEO III and asymmetric B factories at SLAC and KEK should be sufficient for a mixing search. The extra 
path-length due to the Lorentz boost, together with the use of silicon detectors for high resolution position 
measurements, have given the fixed target experiments an advantage over e + e~ experiments (C/a ~ 8— 10). 
The low background around the primary vertex at e + e~ experiments and photoproduction experiments is 
a certain advantage. It is worth pointing out here that at the e + e~ experiments (esp. at asymmetric B 
factory or Z factory) it maybe possible to use B° — > D* + l~v, where the primary (D* + decay) vertex can be 
determined by the l~ together with the slow pion coming from the D* + . In this case, the background level 
around the primary vertex is intrinsically very low. 

However, in the case of D° — ► K + K ~, 7r + 7r~, etc., the requirement on the background level around 
the primary vertex is not so important. In this case, the mixing signature is not a deviation from a perfect 
exponential (again assuming CP conservation), but rather a deviation of the slope from (71 + 72V2. It is 
worth pointing out that there are many advantages with this method. For example, one can use Cabibbo 
favored decay modes, such as D° — ► K~tt + , to measure the average D° decay rate (71 + 72V2. This, 
along with other SCSD CP even (or odd) final states, would allow for valuable cross checks on systematics 
uncertainties. In addition, since we only need to determine the slope here, we do not need to tag the D° nor 
know the primary vertex location. The sensitivity of this method depends on how well we can dertermine the 
slope. Roughly speaking, in the ideal case, the sensitivity to y would be ~ l/\/~N, where N is the number 
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of D a — > K + K~, ir + ir~, etc. events, which means that the sensitivity to mixing caused by the decay rate 
difference (~ y 2 /2) would be close to ~ \/N . For example, a fixed-target experiment capable of producing 
~ 10 8 reconstructed charm events could lower the sensitivity ^] to ~ 10~ 5 — 10~ 6 level for the y 2 term in 
Rmixing = {x 2 +y 2 )/2. 

4.2 Hadronic method B 

In the near future, we should be able to have a good understanding of DCSD in D° — > K + Tr~Tr°, 
D° — > K + ir~Tr + ir~ , etc. modes, then method B will become a feasible way to study mixing and the sensitivity 
should be improved. Just like method A, this method requires very good vertexing capabilities and very 
low background around the primary vertex (this is even more important than in method A, since precise 
knowledge of DCSD is very important here) . In addition, this method requires that the detection efficiency 
(for the mode being searched) across Dalitz plot be quite uniform (at least the detector should have good 
acceptance on the Dalitz plot at locations where DCSD and mixing resonant substructure are different). 
This is necessary so that detailed information on the resonant substructure can be obtained in every corner 
on the Dalitz plot. 

The excellent photon detection capabilities will allow e + e~ experiments to study the D° — > K + ir~ir° 
mode with very low background. From the CLEO II D° — > K + ir~ir° analysis, the detection efficiency across 
the Dalitz plot will have some variations due to cuts needed to reduce background, however, it is still good 
enough to obtain detailed information on the resonant substructure |p5| . Future fixed target experiments 
may have a good chance to study D° — > K + TT~n + Tr~ mode, since the detection efficiency across Dalitz plot 
should be quite flat. The sensitivity that each experiment can reach by using this method depends on many 
things and need to be carefully studied in the future. 

4.3 Hadronic method C 

The sensitivity of this method depends crucially on the particle identification capabilities. Since the D° 
is at rest, the K and ir mesons will have the same momentum, so a doubly misidentified Cabibbo favored 
decay D° — > K~ir + (K~ — > ir~,ir + — > K + ) mimics a D° — > K + tt~ with almost the same D° mass. It is 
worth pointing out here that particle identification is not as crucial to method A as it is to this method (C), 
as far as this particular background is concerned. This is because in method A, the D° is highly boosted, and 
doubly misidentified D° — > K + ir~ decays will have a broad distribution in the D° mass spectrum around the 
D° mass peak; this background can be kinematically rejected with only a small reduction of the efficiency 
for the signal events Q 

Once the sensitivity reaches 0(1O -5 ), one may have to worry about other contributions, such as con- 
tributions from continuum background, contributions from e + e~ — + 27 — > D°D° which may produce C-even 
states where DCSD can contribute [p7f . 



5 Since I came up with this idea the day before the deadline of this paper, the sensitivity has not been 
carefully estimated yet. 

6 It may be possible that good understanding of DCSD can be reached by measuring the pattern of D + 
DCSD decays where the signature is not confused by a mixing component. It is worth pointing out that the 
D + DCSD decays can be studied very well at future fixed target experiments. 

7 The idea is for each K + ir~ candidates, one can invert the kaon and pion assignments and recalculate the 
D° mass, denoted Mm p . If Mfu p lies close to the nominal D° mass, the combination is discarded. This veto 
works as long as the momentum measurement is correct. One can say that excellent tracking capabilities is 
crucial in order to get rid of this background here. 
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4.4 Semi-leptonic method 



The semi-leptonic method usually suffers from large background (except at a r charm factory), the 
traditional method of looking for like sign or fi~fji~ pairs is an example. New ideas are needed in order 

to improve the sensitivity significantly. Some promising techniques have been suggested by Roily Morrison 
and others, and have been discussed in the working group f30| . 

5 Summary 

The search for D a D° mixing carries a large discovery potential for new physics since the D°D° mixing 
rate is expected to be very small in the Standard Model. The past decade has seen significant experimental 
progress in sensitivity (from 20% down to 0.37%). 

In light of the recent CLEO II signal in D° — > K + ir~ , if the mixing rate is close to that of DCSD 
(above 10~ 4 ) , then it might be observed by the year 2000 with either the hadronic or the semi-leptonic 
method, either at fixed target experiments, CLEO III, asymmetric B factories(at SLAC and KEK), or at 
a r-charm factory. If the mixing rate is indeed much smaller than DCSD, then the hadronic method may 
have a better chance as the potentially very small mixing signature could be enhanced by the presence of the 
relatively large "annoying background" DCSD. The design of future experiments should focus on improving 
the vertexing capabilities and reducing the background level around the primary vertex, in order to fully take 
advantage of having the possible DCSD interference. In addition, the very complication due to the possible 
differences between the resonant substructure in many DCSD and mixing decay modes D° — > K + ir~(X) 
could in principle be turned to advantage by providing additional information once the substructure in 
DCSD is understood (method B) and the sensitivity could be improved significantly this way. This means 
that understanding DCSD in D decays could be a very important step on the way to observe mixing. 
Experimenters and theorists should work hard on this. 

In the case of D° — > K + it~(X) and D° — > X + l~, we are only measuring R m i X i ng = (x 2 + y 2 )/2. Since 
many extensions of the Standard Model predict large x — Sm/^/ + , it is very important to measure x and y 
separately. Fortunately, SCSD can provide us information on y. This is due to the fact that decays such as 
D° — > K + K~, 7r + 7r~, occur only through definte CP eigenstate, and this fact can be used to measure the 
decay rate difference y = 7-/7+ = (71 —72)/ (71 +72) alone. Observation of a non-zero y would demonstrate 
mixing caused by the decay rate difference. This, together with the information on R m i x i ng obtained from 
other methods, we can in effect measure x. In this sense, it is best to think of the quest to observe mixing 
(new physics) as a program rather than a single effort. 
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Figure 1: The quest for D°D° mixing. Note that the range in E691 result reflects the possible effects of 
interference between DCSD and mixing, and the CLEO II signal could be due to either mixing or DCSD, or 
a combination of the two. 



Figure 2: The CLEO II signal for D° — ► K + tt~. The D° mass for wrong sign events, (a) for events in the 
AM peak; (b) for events in the AM sidebands. The solid lines are the fits using the corresponding right 
sign mean and a in data. 
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Figure 3: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%. 
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Figure 4: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%, in log scale. 
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Figure 5: The decay time dependence of DCSD and mixing with maximal destructive interference coscf> = 
-1.0. For different a = R mixi ng/RDCSD values: from left to right, a = 100%, 50%, 10%, 5%, 1%,0.5% (with 
Rdcsd = 10~ 2 , this corresponds to R mixing = 10~ 2 , 5 x 10~ 3 , 10~ 3 , 5 x 10~ 4 , 10~ 5 , 5 x 10~ 6 ). 
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Figure 6: The decay time dependence of DCSD and mixing with a = Rmixing/R-DCSD = 10%. For different 
coscj) values: from bottom to top, cos(f> = -1.0,-0.99,-0.96,-0.94,-0.80, -0.6,0.0,0.5,0.7,1.0. The solid 
line is the DCSD term, as a reference line. 
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